Coarse-grained Monte Carlo simulations of the phase transition of Potts model on 

weighted networks 
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Developing effective coarse grained (CG) approach is a promising way for studying dynamics on 
large size networks. In the present work, we have proposed a strength-based CG (s-CG ) method to 
study critical phenomena of the Potts model on weighted complex networks. By merging nodes with 
close strength together, the original network is reduced to a CG-network with much smaller size, on 
which the CG-Ifamiltonian can be well-defined. In particular, we make error analysis and show that 
our strength-based CG approach satisfies the condition of statistical consistency, which demands 
that the equilibrium probability distribution of the CG-model matches that of the microscopic 
counterpart. Extensive numerical simulations are performed on scale-free networks, without or with 
strength-correlation, showing that this s-CG approach works very well in reproducing the phase 
diagrams, fluctuations, and flnite size effects of the microscopic model, while the d-CG approach 
proposed in our recent work [Phys. Rev. E 82, 011107(2010)] does not. 

PACS numbers: 05.50.-|-q, 89.75.Hc, 05.10.-a 



I. INTRODUCTION 

In the last two decades, we have witnessed dramatic 
advances in complex networks research, which has been 
one of the most active topics in statistical physics and 
closely related disciplines [l|-[5[- The central issue in this 
field is to study how the topology of networks influences 
dynamics, such as phase transition, self-organized criti- 
cality and epidemic spreading, etc. Usually, Monte Carlo 
(MC) simulations [GJ have been widely used to study such 
dynamics. However, the sizes of many real-world net- 
works are very large, such as human brain composed of 
about 10^^ neurons and 10^^ synapses [?[, and thereby 
brute-force simulations are quite expensive and some- 
times even become impossible. Phenomenological mod- 
els, such as mean-field description, may capture certain 
properties of the system, but often ignore microscopic 
details and fluctuation effects which may be important 
near some critical points. Therefore, a promising way 
to bridge the gap between the microscopic details and 
system level behaviors is to develop coarse-grained (CG) 
approaches, aiming at significantly reducing the degree 
of freedom while properly preserving the microscopic in- 
formation of interest. 

Recently, several CG approaches have been proposed in 
the literature. Renormalization transformation has been 
used to reduce the size of self-similar networks while pre- 
serving the most relevant topological properties of the 
original ones (8l-[ll|. Gfeller and Rios proposed spectral 
decomposition technique to obtain a CG-network which 
can reproduce the random walk and synchronization dy- 
namics of the original network [l^]- Kevrekidis et al. 
developed equation-free multiscale computational meth- 
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ods to accelerate simulation using a coarse time-stepper 
[l3| , which has been successfully applied to study the 
CG dynamics of oscillator networks 14], gene regulatory 
networks p^ . and adaptive epidemic networks 16]. Nev- 
ertheless, none of the works mentioned above has consid- 
ered critical phenomena in complex networks, which has 
been a frontier topic in the context of network science 

m- 

Very recently, we have proposed a degree-based CG 
[d-CG ) approach to study the critical phenomena of the 
Ising model and the SIS-epidemic model in unweighted 
networks [l^. A local mean field (LMF) scheme was 
introduced to generate the CG network from the micro- 
scopic one. Specifically, we have proposed a so-called 
condition of statistical consistency (CSC) that the CG- 
model should satisfy to guarantee the validity of the CG- 
approach. We showed that the CSC can be exactly ful- 
filled if we merge nodes with the same degree together. 
Extensive numerical simulations showed that our d-CG 
approach does work very well to reproduce the phase 
transition behaviors of the original network, including 
the critical point and the fluctuation properties, but with 
much less computational efforts. Our method also makes 
it feasible to investigate the finite size effects of both 
models, which should be much more expensive and even 
forbidden if we use brute-force methods. However, this 
d-CG approach can only apply to binary networks, i.e., 
each of the link in the network either exists or not, but 
with no weight. As we know, many real-world networks 
are intrinsically weighted, with their links having diverse 
strengths. Examples include the collaboration networks 
fl342l|, airport networks [13, HI], metabolic networks 
24] predator- prey relationship networks [1^, and so on. 
Therefore, a straightforward question is: Can we use CG 
approaches to study the critical phenomena in weighted 
networks? 

To answer this question, in the present work, we have 
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considered the critical phenomena of the Potts model in 
weighted complex networks. The Potts model is related 
to a number of irnportant topics in statistical and math- 
ematical physics [2^ [13] and was successfully applied to 
neural networks, multiclass classification problems, graph 
coloring problem, and so on. It contains a system of cou- 
pled nodes, each of which has p possible states. Only 
when two nodes are in the same state, they have pairwise 
interactions. With the increment of temperature, the 
Potts model undergoes an order-disorder phase transition 
at some critical temperature. For p — 2, Potts model is 
equivalent to the well-known Ising model. Instead of the 
d-CG scheme, we have proposed a strength-based CG (s- 
CG ) approach, where those nodes with similar strength 
are merged together to form a CG-node. Note that in 
weighted networks, it is unpractical to merge nodes with 
exactly the same strength together. By detailed analysis 
of the discrepancy between the Hamiltonian of a CG con- 
figuration and that of its corresponding microscopic con- 
figurations, we show that the s-CG approach can approx- 
imately satisfy the CSC defined on weighted networks. 
Extensive numerical simulations are performed on scale- 
free(SF) networks, without or with strength-correlation, 
showing that our s-CG approach works very well in re- 
producing the phase diagrams, fluctuations, and finite 
size effects of the microscopic model, while the sirnple 
d-CG does not. Compared to our previous work [l8j, 
the present study step forward several important steps. 
First of all, we should note that s-CG is a brand new 
method compared to d-GG and the latter cannot apply 
to weighted networks, although they share some similar 
ideas. Secondly, weighted networks are of more ubiqui- 
tous importance than binary unweighted ones, thus the 
s-CG approach should flnd more applications. What is 
more, we have extended the study from the simple two- 
state Ising model to a more general one, the multi-state 
Potts model. In addition, we have performed error anal- 
ysis in the present study, which clearly demonstrates the 
robustness of our approaches. 



II. COARSE GRAINING PROCEDURE 



lowed by CG-MC simulations to study the dynamic be- 
haviors. The CG-network is simply obtained by node- 
merging, i.e., nodes within the original micro-network 
are merged into a single CG-node C^, where /i = 1, N'^ 
labels the CG-node and N'^ is size of the CG-network. 
Following the LMF scheme used in Ref. [l^, the weight 
of link between two CG nodes /i and u reads. 



for /i = ly, 
for /i =/= u. 



(2) 



The CG-Hamiltonian H can be readily obtained, 
H — Hi + H2 



where 



Hi 
H2 



^ (3a) 



p,u{>p) a 



(3b) 



Herein, Hiifl-i) denote CG interactions inside(among) 
the CG-nodes, respectively. 'r\p,,a stands for the num- 
ber of a-state micro- nodes inside Cn- Since there are 



possible distinct pairs of a-state micro-nodes 



inside , and each pair has a weighted coupling w^^ , the 
CG-interactions among all the a-state nodes inside is 
given by 



(a) 



VpAVp^a - 1) 



Summation this over all CG-nodes fi and states a gives 
the result in Eq.(3a). Eq.(3b) can be interpreted in a 
similar way. Note that Eq.(3) are closed at the CG level, 
i.e., as long as one has constructed the CG-network, w 
and H are then both well defined, based on which one 
can perform CG-MC simulations without going back to 
the micro-level. 



CG Potts Model 



B. CSC: Condition of statistical consistency 



In this paper, we consider the p-states Potts model on 
a weighted network consisted of TV nodes, whose Hamil- 
tonian is given by 



(1) 



wher 



is the weight on the edge connecting a pair 



of nodes i and j {wij = if the nodes i and j are not 
connected). ai{— 1, • • • ,p) denotes the state of node i, 
^ai,aj = 1 if a,; — Oj and otherwise. 

To setup the CG-Potts model, one needs to obtain 
the CG-Hamiltonian defined on the CG-network, fol- 



The above procedure tells us how to calculate the CG- 
Hamiltonian if we already have the CG-network. How- 
ever, which Qfj^ nodes are merged together to form a 
CG-node is yet not determined. Generally speak- 
ing, one may construct the CG-network deliberately, for 
instance, one may simply generate N'^ values, obey- 
ing X]^=i = and then just randomly merge 
micro- nodes to form C^. Therefore, an important ques- 
tion arises: How to guarantee that the CG-model can re- 
produce the dynamics of the corresponding microscopic 
model correctly? 

We address this problem by extending the so-called 
CSC as proposed in [l^. We introduce ff^ = 
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to denote the state of and fj = 

{'iffj.}fj.=i N'= to denote the configuration of the CG- 

network. Note that a given CG configuration i] corre- 
sponds to many microscopic configurations, which defines 
the degeneracy factor g(jj). In the equilibrium state of 
the CG-model, the probabihty of finding a given CG- 
configuration fj is given by the canonical distribution, 
i.e., 



Pcoiv) = gir])e 



-H/ksT 



/Z, 



where Z = ^ffPcoiv) is the CG partition function. It is 
important to note, however, that ^(17) can be calculated 
exactly from the equilibrium distribution of the micro- 
model. 



Pn 



,-H/kBT 



/z, 



where Z is the partition function of the micro-model, and 
the prime means summation over all the microscopic con- 
figurations that contribute to rj. Since we are interested 
in the equilibrium phase transition behavior of the Potts 
model, we thus assert that for the CG-model to be sta- 
tistically consistent with the micro- model, pcciv) ^^'^ 
Pmicro('i7) must bc cqual, i.e., the CSC reads 



'H/kBT 



/z. 



(4) 



C. s-CG Scheme and error analysis 



2 (g*^ + 5si){Si_, + 5sj) 



SN 



— E 



{Sf, + Sst){Si, + 6sj) 



SN 
SN 



(6a) 



(6b) 



Herein, we have written Si = 5*^ -|- Ssi, with S^ ~ 
~ Siec being the mean strength within C^. il^ = 



9m 

^here (^s^)^ ^ ^ E^ec,, is the variance 

of strength within C^. In the first equation, we have 
used the fact that (Eiec^ '^■^O^ = 2 Ei<jec,, '^^i'^^J + 
Eigc^C'^'Sj)^ = 0- The second equation holds simply be- 
cause Ei6C^,i6C„ -^^f^^J = (Eisc^ '^Si)(EjGC„ -^Si) = 0- 
Substituting Eq.® into Eq.Q, we can get 



^1 - -^E^^(i-".)E 



Vti,a{Vt^,a - 1) 



^2 = ^ E 



(7a) 
(7b) 



In the present work, we propose a s-CG scheme to 
construct the CG- network, i.e., nodes with same or sim- 
ilar strengths are merged together to form a CG-node, 
where the strength Si of node i is defined as Si — Wij 
[23 , [28|. In the following, we will show that if nodes 
inside each CG-node have same strengths, the CSC will 
hold exactly within the ANA. In addition, if the strengths 
within are nearly the same, the CSC can also hold ap- 
proximately. 

In the literature, ANA [l^,!!!-!!! has been widely used 
to study the ensemble averaged dynamics of complex net- 
works and proved to be successful. ANA assumes that 
one can replace the dynamics on a given network by that 
on a weighted fully connected graph with connectivity 
Aij — didj/{DN), where di (dj) denotes the degree of 
node i (j) and D is the mean degree of the network. 
Analogously, in weighted networks link weight can be 
expressed as 

Wij = SiSj/{SN) (5) 

where S is the mean strength of the network. Substi- 
tuting Eq.lO into Eq.®, the adjacency matrix of the 
CG-network now reads. 



To compare the CG-Hamiltonian with the microscopic 
one, we now group the micro- nodes with same state a 
inside as C^,q.. Clearly, the size of C^,q is ilp,,a- As in 
Eq.Q, we can also split the micro-Hamiltonian H into 
two parts , 

H = Hi+H2 (8) 

where Hi and denote energy contributions from intra 
and inter the CG-nodes respectively. With ANA, and 
noting the fact only nodes with same states have inter- 
actions at the micro-level, one has 

Hi = -SN 2 

^2 = - E E E (9b) 



Following similar steps to obtain Eq.Q, we may also 
write Si = S'^^q + Ssi (here node i belongs to the group 
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C^,q) and Eqs.® change to 



EE E 

VlJ.Ar]l^,a - 1) 



SN 



such that Sfj.,a averaged over time is close to 5*^. Hence 
{H — H)/H averaged over long time could be small. Note 
that if we merge nodes randomly, ^ 1 and ^i^^Q «C 1 
will be violated and the above reasoning should fail. We 
thus conclude that the pratical s-CG approach, by merg- 
ing nodes with similar strength together, can satisfy the 
CSC approximately. 



(10a) 



— E E E 



SN 



SN 



(lOb) 



/i,i/(>/j) a 



Here = sCivlT-^) '^^^'^^ (^^')p,a = 

Siec i^Si)"^ is the variance of strength within the 
group of nodes C^^q. 

Comparing Eq.Q with Eq. (fTO|) . the discrepancy be- 
tween the CG-Hamiltonian and the micro-Hamiltonian, 
is given by 



Hi -Hi 



SN 



III. NUMERICAL RESULTS 

To show the validity of our s-CG approach, we per- 
form extensive simulations on weighted SF networks. SF 
networks are much heterogeneous and serve as better can- 
didates to test our method than other homogeneous net- 
works, such as small-world or random networks (other 
types of complex networks have also been investigated, 
the qualitative results are the same and not shown here) . 
We first generate a regular (unweighted) SF network by 
using the Barabasi-Albert (BA) model [s^ with power- 
law degree distribution P(fc) ~ k~^. To convert this un- 
weighted SF network into a weighted one, we use the 
algorithm as proposed in Ref. [33: The weight of a 
link between node i and j (1 < i,j < N) is given by 
Wij = {jj + /2, where 6* is a tunable parameter. Note 
that 9 = corresponds to an unweighed network. 

The MC simulation at the microscopic level follows 
standard Metropolis dynamics: At each step, a micro- 
node is randomly selected and its state is randomly up- 
dated with an acceptance probability min(l, Q-^H/ksT-^^ 
where A_ff is the associated change of the micro- 
Hamiltonian, ks is the Boltzmann constant, T is the 
temperature. In the present work, we set fcs = 1. Sim- 
\/{jj,,a) and S*^ = S^,a, hence Hi = Hi ilarly, during each CG-MC step, a CG-node is ran- 
domly chosen with probability proportional to its size g^. 
The probability for the process that an a-node changes 
to a /3-node, with correspondingly rj^^a ^ Vfj..a ~ 1 and 
Vfj.,^ Vti,i3 + ^, is given by ?7^^„ min(l, e~^^/''«^), where 
AH is the change of CG-Hamiltonian during this pro- 
cess. Since N'^ can be much smaller than N, the CG-MC 
is expected to be much faster and memory-saving than 
the micro-level MC simulation. 

The collective state of the system is described by 
the total magnetic moment M = J2 



fi a 

[sl{i-n^)-slM-^^.A] (iia) 
H2-H2 = -2^2^ ^ ^ib) 



Obviously, for the exact s-CG algorithm where all 
the nodes inside a given CG-node have same strengths. 



and H2 = H2- In this case, all those microscopic config- 
urations contributing to a CG-configuration 77 have ex- 
actly the same Hamiltonian H , which also equals to the 
CG-Hamiltonian H. Since the constrained summation 
contains exactly 5(17) items, the numerators on both 
sides of Eq.dH) are exactly equal, i.e., g{fj)e~^ /^'^'^ = 
e"^/*"'^"^. Since we can also write the microscopic 
partition function as Z = ^f^i^' it is readily 

to show that the two partition functions equal, Z = Z . 
Therefore, the CSC, Eq.(|l]), exactly holds. 

However, we should note that for a weighted net- 
work, the exact s-CG method is not practical, since 
the strength of a given node is generally not an integer. 
Therefore, usually one can only merge nodes with close 
strength together. Let us analyze Eq. pTj) again. The 

factor Vt^ scales as g-i^'^ , hence if we merge many nodes 
with similar strengths together, fi^ ^ 1 is expected to be 
true. One may also expect that il^^a 1 for the same 
reason. Therefore, the discrepancy between H and H 
mainly depends on the difference between S*^ and 6*^,0. 
Here, we note that the nodes with a-state fiip with time 
during the simulation. In the equilibrium state, one ex- 
pects that C^,a may scan throughout for many times, 



Mf,,a\, where 



p-i 



(^ = l,---,iV'^) denotes the 



component of the magnetic moment within C^. With 
increasing temperature T, the Potts model undergoes a 
phase transition at some critial temperature Tc from an 
ordered state, where M ~ 0(1) is strictly nonzero, to a 
disordered state with M ~ 0. We use the similar s-CG 
approach to construct the CG-network with different N'^ 
and compare the results obtained from CG-MC simula- 
tions with those of micro-MC simulations. 

To begin, we show the results in Fig.l for 61 = 0, where 
the network are essentially unweighted and the s-CG ap- 
proach is identical to the d-CG . Fig. 1(a) and 1(b) show 
the moment M and susceptibility x = I3N{{M'^) - (M)^) 
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as functions of T, respectively. The susceptibility is 
related to the variance of the total magnetization ac- 
cording to the fluctuation-dissipation theorem. Appar- 
ently, our results (empty squares and solid circles) are in 
excellent agreements with the micro-level counterparts 
(solid lines). As comparisons, we have also shown the 
results obtained by a random-merging (RM) CG-model 
(dotted lines) and the heterogeneous mean field theo- 
ries (HMFT)[3l| (empty triangles). Here, the RM model 
means that one simply merge N/N'^ randomly selected 
nodes to form a CG-node. Evidently this random scheme 
fails to reproduce the microscopic behaviors at all. The 
results of the HMFT are obtained by numerically solv- 
ing the self-consistent equations of order parameter js^. 
We find that the HMFT can predict the curve of M - T 
quite well, however, it fails to predict the curve of x ~ T. 
Strikingly, even when the original network is reduced to 
one with only 16 CG-nodes, the CG model still faith- 
fully reproduces the phase transition curves and fluc- 
tuation properties. Since N'^ is largely reduced com- 
pared to N, a considerable speed-up of CPU time can 
be achieved which makes it feasible to study system size 
effects. Fig. 1(c) plots Tc as a function of In A^, obtained 
by our GG method with N'^ = 64. Tc is determined 
as the location of the peak in the x ^ T curve, see 
Fig. (lb). The dependence is linear with a slope ~ 1.68, 
which agrees rather well with a theoretical prediction 
Tc/lnN = ^ =i 1.67 0], where S is the average node 
strength in the network. 




FIG. 1: (color online), (a)-(b) M and x as functions of T (in 
unit of J/ks) for the Potts model on unweighted SF networks 
(9 — 0), obtained from brute- force MC simulation (solid line), 
HMFT (triangle), random- merging CG(dotted line), and the 
s-CG (square and circle). A'' — 16384, p = 3. (c) Dependence 
of Tc on the network size TV obtained by the s-CG approach 
with fixed A'''^ — 64. All the networks have fixed mean degree 
D = 20. The error bars (not shown) are smaller than the 
symbol sizes. 

For 6* 7^ 0, the networks are weighted. Here we take 
6 ~ 2.4 as an example to ensure the heterogeneity of 
the link weights. Figure(2a) and (2b) show M and x 
as functions of T respectively. As in Fig. 1(b), the peak 
in X locates the critical point Tc- Clearly, the s-CG re- 
sults (solid circle) are still in excellent agreements with 
the MC results (solid lines), however, the c?-CG (solid 
squares) [11] and RM-CG (dotted lines) both fails. For 
such weighted networks, the dynamic equations of HMFT 
is not available either. Thus, for such weighted net- 
works, our s-CG approach is the only promising CG ap- 
proach so far. In Fig. 2(c), we have also shown the de- 



pendence of Tc on the network size. Apparently, there 
is also a linear dependence between Tc and InA^ with 
the slope being about 1.288. As mentioned in the last 
paragraph, this slope depends on the average strength S. 
For a weighted network, one may estimate S by {wij)D, 
where (wij) ~ J^x^/Adx — j(jippj'2^^^ ■ Substituting 
D — 20, 6 — 2.4 and p = 3 to these formula, we ob- 
tain Tc/\nN = ^ ~ 1.293, which is consistent with the 
simulation value. 




T T InN 



FIG. 2: (color online), (a)-(b) M and x as functions of T 
for the ferromagnetic Potts model on a weighted scale-free 
network with mean degree D = 20. N ^ 1024, p = 3, 6* = 2.4 
and A^'^ — 16. (c) Dependence of Tc on the network size A*'. 
The error bars are omitted for clarify since they are smaller 
than the symbol sizes. 

In real- world networks, correlation is an ubiquitous 
feature. For instance, social networks show that nodes 
with large degrees tend to connect together, a prop- 
erty referred to as "assortative mixing" [35i] . In contrast, 
many technological and biological networks show "disas- 
sortative mixing", i.e., connections between high-degree 
and low-degree nodes are more probable [36l . [37j . Previ- 
ous studies showed that correlations may play important 
roles in network dynamics ^35l-l39j. In the present work, 
we have used our s-CG method to study the phase tran- 
sition of Potts model on correlated networks, which can 
not be studied by the HMFT which assumes no degree 
correlation. To characterize the assortative property of 
the weighted network, a strength correlation coefficient r, 
an extension of the degree correlation f35j , can be defined 
as 

r=((s,s,)-(s,)(s,))/((s2)-(s,)2). (12) 

Here Si and Sj are the strengths of the two end-nodes 
of an edge, r is zero for networks with no strength- 
correlation, such as BA-SF networks, and positive or neg- 
ative for assortative or disassortative mixing networks, 
respectively. 

Figure 3(a) shows Tc as a function of r, obtained from 
our s-CG approach and micro-MC simulations for = 0. 
Again, the fits between CG-MC and MC are good. Fig- 
ure 3(b) shows the effects of correlated network size on Tc- 
Interestingly, we find that the linear dependence between 
Tc and InA^ is lost for correlated networks. For assor- 
tative(disassortative) networks Tc grows monotonically 
much faster(slower) than IniV, respectively. In other 
words, the ordered state in an assortative(disassortative) 
network is harder(easier) to be destroyed with increas- 
ing temperature than in an un-correlated network. This 
is understandable since a 'hub'-node in the network is 
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FIG. 3: (color online). Phase transition behaviors of the Potts 
model on unweighted correlated networks, (a) Tc plotted as a 
function of the network correlation coefficient r, obtained via 
CG-MC and MC simulations. iV = f024 and A^"^ = 64. (b) 
Dependence of Tc on the network size TV. All the networks 
have fixed mean degree D = 20. 

more difficult to change its state than a 'leaf-node due 
to larger energy barrier. In an assortative network, hub- 
nodes are connected together, such that they tend to 
freeze into a local ordered state which is stable to ther- 
mal fluctuations. For a disassortative network, a hub- 
node is usually connected to many leaf-nodes. Since leaf- 
nodes can change state easily, the 'alone' hub-node is 
more likely to change state with the help of their 'boiling' 
neighbors. Therefore, assortative correlations tend to in- 
crease Tc as observed here. 




5 10 15 2 4 6 8 10 12 14 



T T 

FIG. 4: Magnetization M and susceptibility x ^ functions 
of temperature T for the Potts model on weighted correlated 
networks. Symbols and lines correspond to the CG-MC and 
micro-MC simulation results, respectively. Other parameters 
are same as in Fig. 2. 



In Fig. (4), the magnetization M and susceptibility x of 
the ferromagnetic Potts model on weighted networks are 
plotted as functions of temperature T at different correla- 
tion coefficient r, obtained from our s-CG approach and 
micro-MC simulations. Again, the agreements between 
CG-MC and MC are excellent, further demonstrating the 
validity of our method. 



IV. CONCLUSIONS 

In summary, we have developed a stength-based s-CG 
approach for coarse-graining study of the phase transi- 
tion of the Potts model on weighted networks. We have 
utilized a mean-field scheme to generate the connectiv- 
ity of the CG-network and derived the CG-Hamiltonian. 
To address the problem how to guarantee the validity 
of the GG- model, we have proposed the so-called CSC, 
which requires that the probability to find a given CG- 
configuration in the equilibrium state, calculated from 
the CG-model, should be the same as that calculated 
from the original microscopic model. We show, by per- 
forming error analysis, that our s-CG approach, by merg- 
ing nodes with close strengths together, holds the CSC 
approximately with ANA. Detailed numerical simula- 
tions demonstrate clearly that our s-CG approach can 
reproduce the microscopic MC simulation results very 
well, not only for the onset of phase transition, but also 
for the fluctuations and system size effects. 
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